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Based on the unified Hauser—Feshbach and exciton model, which can describe the particle emission processes 
between discrete energy levels with energy, angular momentum, and parity conservations, a statistical theory of 
light nucleus reaction (STLN) is developed to calculate the double-differential cross-sections of the outgoing 
neutron and light charged particles for the proton-induced °Li reaction. A significant difference is observed 
between the p + °Li and p + “Li reactions owing to the discrepancies in the energy-level structures of the targets. 
The reaction channels, including sequential and simultaneous emission processes, are analyzed in detail. Taking 
the double-differential cross-sections of the outgoing proton as an example, the influence of contaminations 
(such as 'H, “Li, '*C, and 180) on the target is identified in terms of the kinetic energy of the first emitted 
particles. The optical potential parameters of the proton are obtained by fitting the elastic scattering differential 
cross-sections. The calculated total double-differential cross-sections of the outgoing proton and deuteron at 
Ep = 14 MeV agree well with the experimental data for different outgoing angles. Simultaneously, the mixed 
double differential cross-sections of He and a are in good agreement with the measurements. The agreement 
between the measured data and calculated results indicates that the two-body and three-body breakup reactions 
need to be considered, and the pre-equilibrium reaction mechanism dominates the reaction processes. Based on 
the STLN model, a PLUNF code for the p + Li reaction is developed to obtain an ENDF-6-formatted file of 
the double-differential cross-sections of the nucleon and light composite charged particles. 
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I. INTRODUCTION 


Nuclear data (including the cross-sections of all types of 
reaction channels, differential cross-sections, and double- 
differential cross-sections) of nucleon-induced ®:7Li reactions 
are important for various applications, such as the compact 
accelerator-driven neutron source and the International Fu- 
sion Materials Irradiation Facility (IFMIF), as introduced 
in Refs. [1-7]. A lithium-glass scintillator with natural 
lithium is used as a small-angle neutron scattering spectrome- 
ter (SANS) for material research at the China Spallation Neu- 
tron Source (CSNS) to satisfy both the neutron detection effi- 
ciency and gamma elimination requirements [8]. “Li enrich- 
ment and °Li inventory are important influencing factors in 
molten salt fast reactors (MSFRs) [9]. In recent years, lithium 
has received considerable attention in astrophysics because it 
can lead to a better understanding of the evolution and forma- 
tion of the universe. The standard big-bang nucleosynthesis 
(BBN) theory overestimates the primordial “Li abundance by 
a factor of approximately three to four, that is, the cosmolog- 
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ical lithium problem [10]. The reaction °Li(p, y) “Be plays 
an important role in the consumption of °Li and the forma- 
tion of "Be, and Be will eventually decay into “Li at the 
end of big-bang nucleosynthesis [11]. In an atmosphere of 
approximately 1% of giant stars, there exists an anomalous 
elevation in lithium abundance, directly contradicting the ex- 
pectations calculated by conventional stellar evolution mod- 
els. An extremely Li-rich giant (possibly newly enriched) 
and a rigorous investigation of its evolutionary stage are def- 
initely important [12]. Furthermore, an in-depth study of the 
nucleon- ° "Li reaction will enhance our understanding of 1p- 
shell light nuclear reactions[ 13-15]. 


Although ĉLi and ‘Li differ by only one neutron, they ex- 
hibit significant differences. The natural abundance of Li is 
only 7.5%, whereas that of “Li is as high as 92.5%. The en- 
ergy levels of “Li and “Li are different both in terms of their 
energy values and widths. Additionally, every energy level 
of the °Li nucleus has a spin with an even multiple of 1/2, 
whereas the spin of the “Li nucleus is an odd multiple of 1/2 
[16]. Li is easily captured by neutrons, producing "Li and 
high-energy gamma rays, whereas Li is more likely to un- 
dergo a decay. Moreover, the capture cross-section of °Li is 
smaller than that of “Li in the thermal neutron energy range, 
making °Li a common material for slow neutron shielding in 
certain nuclear applications. Differences in the energy level 
diagrams of Li and ‘Li can lead to variations in their de- 
cay properties and the particles emitted from the individual 


levels [1-6]. Therefore, the results of the n+®”Li reactions 
are significantly different, as shown both experimentally and 
theoretically [1, 2, 7]. In our previous study, we successfully 
predicted the double-differential cross-sections of the charged 
particles for the p+"Li reaction using the statistical theory of 
light nucleus reaction (STLN) model [15]. On this basis, 
an ENDF-6-formatted file of the double-differential cross- 
section was obtained. File-6 (file of the double-differential 
cross-section), one of the important files of the nuclear reac- 
tion database, is recommended when the energy and angular 
distributions of the emitted particles must be coupled, when 
it is important to provide a concurrent description of neutron 
scattering and particle emission, when so many reaction chan- 
nels are open that it is difficult to provide separate reactions, 
or when accurate distributions of the charged particle or resid- 
ual nucleus are required for particle transport, heat deposition, 
or radiation damage calculations [17]. However, there are 
still no publications describing the double-differential cross- 
sections of outgoing particles for the p+°Li reaction. 


Quantitatively describing all the physical quantities is chal- 
lenging because of the limited availability of comprehensive 
experimental data for the p+°Li reaction. For the p+®°Li re- 
action, there are only a small number of experimental partial 
cross-sections of the (p, y), (p, el), and (p, 3He) channels 
[18-20]. Furthermore, numerous elastic scattering angular 
distributions have been measured at various incident energies 
[20, 21]. Fortunately, the double-differential cross-sections 
of the outgoing proton, deuteron, and mixture of “He and 
a for the p+°Li reaction were measured in 1989 and 1991 
[22, 23]. This makes it possible to validate the theoretical 
calculations. The evaluated partial cross-sections of (p, el), 
(p,° He) and (p, x) are available in nuclear reaction databases 
such as ENDF/B-VIII.0 [24] and JENDL-5 [25]. However, 
the double-differential cross-sections of the outgoing neutron 
and charged particles were not included in these databases be- 
cause the effects of the secondary particle emission processes 
between the discrete levels and cluster separations in light nu- 
clear reactions have not been considered. In ENDF/B-VIII.0, 
the cross-sections of (p, el), (p, 3He), and (p, x) were de- 
rived from the R-matrix analysis [24, 26], while in JENDL-5, 
multichannel R-matrix fitting was used to evaluate the exper- 
imental data in the incident proton energy range from 10 keV 
to 3 MeV. The CCONE code [25] was employed to calcu- 
late the differential cross-sections of the emitted particles in 
the incident proton energy range from 3 to 200 MeV. Subse- 
quently, proton-induced reactions on ® "Li were also calcu- 
lated using the CDCC method in 2013 [5], but the double- 
differential cross-sections of the outgoing proton, deuteron, 
and triton were not provided. A possible reason for this is that 
the sequential secondary particle emission processes were not 
considered in the CDCC model [6]. 


This study provides an in-depth analysis of the p+°Li re- 
action based on the STLN model with energy, angular mo- 
mentum, and parity conservations. The Coulomb barriers of 
both the incident and exit channels are considered for differ- 
ent charged particles. Taking an impure target as an example, 
the influence of contamination is described in terms of the 
kinetic energy. The double-differential cross-sections of the 


total outgoing proton, deuteron, and the mixture of ?He and a 
are calculated by analyzing the reaction channels at different 
incident energies. The calculated results at Ep = 14 MeV 
agree well with the available experimental data. 

Section II briefly introduces the theoretical model used in 
this study. The reaction channels of the p + °Li reaction be- 
low 20 MeV are analyzed in detail in Sect. III. Section IV 
provides a comparison between the calculated results and the 
experimental data, along with the corresponding analysis. A 
summary is provided in the final section. 


Fig. 1. Schematic of the momentum conservation in the nuclear re- 
actions in LS. 
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Fig. 2. (Color online) Measured and calculated total double- 
differential cross-sections of the outgoing proton for a proton- 
induced impure °Li reaction with an angle of 60° at Ey = 14 MeV 
in LS (a). The experimental data are obtained from Ref. [23]. The 
black narrow bands of the below panel (b) represent the contribu- 
tions of the outgoing proton from the discrete levels (only illustrated 
from ground state to the fifth energy level as marked) of targets °Li 
and contaminants (including 1H, “Li, RC, and 160), 


Il. THEORETICAL DESCRIPTIONS 
A. Theoretical Frame 


Based on the unified Hauser—Feshbach and exciton model 
[27], which can describe the particle emission processes be- 
tween the discrete energy levels with energy, angular momen- 
tum, and parity conservations, a statistical theory of light nu- 
cleus reaction (STLN) is developed to describe the mecha- 
nism of the nucleon-induced light nucleus reaction [13, 14]. 
A considerable amount of experimental data, with a focus 
on double-differential cross-sections (such as the neutron- 
induced reactions on Li [1], “Li [2], Be [28, 29], '°B [30], 
11B [31], ‘°C (27, 32, 33], 14N [34], 160 (35, 36], and °F 
[37] as well as the proton-induced reactions on “Li [15] and 
IBe [14, 38]), has been reproduced very well. 

The cross-sections of the first emitted particles from the 
compound nucleus to the discrete energy levels of the first 
residual nuclei can be expressed as 
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where PIT (n) is the occupation probability of the n-th ex- 
citon state in the jm channel, and (j and m denote the an- 
gular momentum and parity in the final state, respectively). 
PIT (n) can be obtained by solving the j-dependent exciton 
master equation under the conservation of angular momen- 
tum in pre-equilibrium reaction processes [39]. QÍ" (n) is 
the occupation probability of the equilibrium state in the jm 
channel. wer p, (n, E* , Ef, ) is the emission rate of the first 
emitted particle mı in thie n-th exciton state with an outgo- 
ing kinetic energy £, in the center-of-mass system (CMS), 
and WÅ” (n, E*) is the total emission rate in the n-th exciton 
state. wie ky (E*, Em, ) is the emission rate of the first emit- 
ted particle my in the equilibrium state with the outgoing ki- 
netic energy €$, in the CMS, and W” (E*) is the total emis- 
sion rate in the equilibrium state. E* is the excited energy of 
the compound nucleus, and gÍ” ( E) is the absorption cross- 
section in the jm channel. The first term of Eq. (1) inside 
the brackets represents the pre-equilibrium process, which is 
the predominant process in | p-shell light nuclei reactions in- 
duced by the nucleon. The second term of Eq. (1) inside the 
brackets describes the equilibrium process. 

The cross-section of the second outgoing particle from the 
discrete energy levels of the first residual nucleus to the dis- 
crete energy levels of the second residual nucleus can be ex- 
pressed as 


Oki ke (n, mı, Mmo) = Ok, (n, mı) i Rii rra (Ex,), (2) 


where op, (n, mı) is the cross-section of the first emitted 
particle, m; expressed in Eq. (1), and Rf1>**2(E,.,) is the 


branching ratio of the second outgoing particle mz from the 
energy level Ep, of the first residual nucleus M; to the energy 
level Ek, of the second residual nucleus M3. 

For the simultaneous emission process, for example, a + 
A — bı + bg + bz (also named after the breakup reaction), 
the total kinetic energy Ec in the CMS can be expressed as 


Ma 
E =A 
CS ap +Q, (3) 


where ma and M4 are the masses of projectile particle a and 
target A, respectively. Fy, is the incident energy, and Q is 
the reaction energy. In terms of the Ohlsen theory [40], the 
momentum distribution function of the bı particle with mo- 
mentum k’; in the CMS is expressed as 


(4) 


12 
Where m; denotes the mass of the b; particle, and c; = — 


Wis 
denotes the kinetic energy of the b; particle in the CMS. The 5 
functions represent momentum conservation and energy con- 
servation. Assuming that the momentum distribution function 
p is a constant, that is, uniformly distributed in the momen- 
tum space. A momentum transformation can be performed 
as 
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one can obtain the following expressions: 
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After performing double-vector integrations, Eq. (7) can be 
expressed as 
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In terms of the isotropic energy spectra in the CMS, as 
mentioned earlier, the energy spectra of bı particle can be ex- 
pressed using the momentum distribution function, i.e., 


an | Ni(ex)des = f hëna, 
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Thus, the energy spectra of the bı particle with outgoing ki- 

netic energy £1 or momentum k’, in the CMS is expressed 

as 


Ny (£1) = C34/ £1 (eT * — €1), (10) 
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maximum outgoing kinetic energy of the bı particle. C3 = 
z is the normalized constant given by f Nı(e1)de1 = 
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1. 
Based on the isotropic energy spectra assumption in the 
CMS for the three-body breakup reaction, that is, the spec- 
tra for all azimuth angles are identical, the double-differential 


cross-section of bı particle in the CMS can be expressed as 
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Similarly, the double differential cross-sections of bz and b3 
can also be expressed in the same form. 

Ec and dQc = dcos cdyc are used to represent the ki- 
netic energy and azimuth angle of the b; particle in the CMS, 
respectively, without losing generality. Similarly, € and 
dQ, = dcos 6dr, are used to represent the kinetic energy 
and azimuth angle of the b; particle in the LS, respectively. 
There is yc = yz both in the measurements and theoretical 
calculations. The following relationship universally exists for 
the double differential cross-section in the CMS and LS: 
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According to the Jacobian determinant of the coordinate 
transformation, the double-differential cross-section of the b; 
particle in the laboratory system (LS) is expressed as [41] 
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Eqs. (10-13) have been validated in our previous studies, 
such as n+°Li >” Li* > n+d +a [1], nt7Li 38 Li* > 
n +t +a [2], and p+’Li >® Be* > p + t +a [15]. 


B. Coulomb Barrier 


Owing to the effect of the Coulomb barrier [42, 43], the 
kinetic energy of the first outgoing charged particle <7,,, must 
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Fig. 3. (Color online) Total double-differential cross-sections 
of the outgoing proton for the p + °Li reaction with angles of 
20°, 30°, 40°, 50°, 60°, and 70° at Ep = 14 MeV in LS. The black 
points denote the experimental data taken from Ref. [23]. The red 
solid lines denote the calculated results. The abrupt peaks represent 
the contributions from the contaminants, such as +H, “Li, Ro and 
160, The different outgoing angles are indicated in the figure. 


be higher than that of the Coulomb barrier Vcou, that is, 
Em, > Vcoul- The reduced penetration factor calculated using 
the optical model potential must be 0 if ef, < Vcou. As- 
suming a spherical nucleus [13], Vcouw can be approximated 


as 
e ZM, Z, 
Mı mı 


a eee ae (14) 
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where Zm, Am, and Zm, Ám, are the charge and mass 
number of the residual nucleus and first outgoing charged par- 
ticle, respectively. rc=1.2—1.5 fm is the charge radius. For the 
charge radii of the proton, deuteron, triton, 3He, a, and °He, 
the experimental data presented in Ref. [44] are used. 

Therefore, the incident energy Ep must satisfy the follow- 
ing equation for open reaction channels: 


Mc 
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where Mc, Mr, and M; are the masses of the compound, tar- 
get, and first residual nuclei, respectively, after the first par- 
ticle is emitted. Æp, denotes the excited energy of the k-th 
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Fig. 4. (Color online) Same as Fig. 3, but at outgoing angles of 
80°, 90°, 100°, 110°, 120°, and 130°. 
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Fig. 5. (Color online)Same as Fig. 3, but at outgoing angles of 
140°, 150°, 160°, and 165°. 


discrete level of the first residual nucleus. By is the binding 
energy of the first particle emitted into the compound nucleus. 
By is the binding energy of the incident particle in the com- 
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Fig. 6. (Color online) Partial double-differential cross-sections of 
the outgoing proton from the p+°Li reaction with an outgoing angle 
of 60° at Ep = 14 MeV in LS. The black points denote the experi- 
mental data taken from Ref. [23], and the red solid line denotes the 
calculated total double-differential cross-sections. The yellow dash- 
dotted line denotes the partial spectra of the emitted proton from the 
direct three-body breakup through “Be* —> p+d+qa. The pink dash- 
dotted lines denote the partial spectra of the first outgoing proton 
from the compound nucleus to the fifth excited energy levels of the 
first residual nucleus, °Li*. The green dash-dotted line denotes the 
contribution of the reaction channel (p, np)’Li* — (p, n2p)a from 
the first excited energy level of °Be* to the first excited energy level 
of °Li*, which can break up into d+a. The brown dash-dotted lines 
denote the partial spectra of the second outgoing proton from the 
fourth and fifth excited energy levels of °Li* to the ground state of 
°He*. The purple dash-dotted lines denote the partial spectra of the 
second outgoing proton from the ground state and the first excited 
energy level of °Li* to the ground state of “He. 


pound nucleus. Clearly, the Coulomb barrier can significantly 
affect the open reaction channels. 


C. Particle Identification 


To conveniently describe the expression of the A(a, b) B re- 
action, some quantities are defined as follows: A, B, a, and b 
are the target nucleus, residual nucleus, incident particle, and 
emitted particle, respectively. The target nucleus A is fixed 
in the LS, so the kinetic energy and momentum are 0, respec- 
tively. The reaction value Q is thus expressed as follows: 


Q = Ey + Ep — Ea, (16) 
where F», Ep, and Ea are the kinetic energies of b, B, and a 
in the LS, respectively. A schematic of momentum conserva- 


tion in the nuclear reactions in the LS is shown in Fig. 1, and 
the momentum conservation is expressed as 


Pa = PB + Po, (17) 
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Fig. 7. (Color online) Same as Fig. 3, but for outgoing 7He and a 
particles. The experimental data are taken from Ref. [22]. 


where Pa = V2Ma Ea, Pb = V2m Ep, and PB=VvV 2mpgEpg 
are the momenta of a, b, and B in the LS, respectively. Fur- 
thermore, Ma, Mp, and mp are the masses of a, b and B, 
respectively. 


From Eqs. (16) and (17), the kinetic energy of the first 
emitted particle b in the LS can be expressed as 
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I. ANALYSIS OF THE REACTION CHANNELS 


For the proton-induced 6Li reaction, reaction channels the- 
oretically exist at an incident energy E, < 20 MeV in terms 
of the reaction threshold energy Ep as follows: 


(p, y) Be, Q = +5.606MeV, Em = 0.000MeV 
(p, n)Be, Q = —5.071MeV, Ep = 5.9206MeV 
(p, pÊLi, Q= 0.000MeV, En = 0.000MeV 
(p,? He)a, Q=+4.019MeV, Em = 0.000MeV 

i (p,d)°Li, Q = —3.442MeV, En = 4.0187MeV 

p+°Li >” Be* > (p, npřLi, Q = —5.666MeV, Em = 6.6153MeV a) 

(p, pn)řLi, Q=—5.666MeV, En = 6.6153MeV 
(p, 2p)5He, Q = —4.594MeV, En = 5.3637MeV 
(p, pd)'He, Q = —1.475MeV, En = 1.7221MeV 
(p,dp)'He, Q = —1.475MeV, En = 1.7221MeV. 


From Eq. (19), one can see that there are obvious differ- 
ences in the p+" Li reaction, as shown in Eq. (20) in Ref. [15]. 
For the first particle emission processes, the p+°Li reaction 
lacks the (p, t) and (p, a) channels. For the second-particle 
emission process, the p+°Li reaction lacks the (p, n°He), (p, 
pt), (n, 2p), (n, 2d), and (n, npd) channels. Furthermore, the 
threshold energies of the same channels are different. 


Considering the energy, angular momentum, and parity 
conservations in the particle emission processes, the reaction 
channels for the first particle emission are as follows: 


n +ê Be* (kı = gs, 1), 
p+®Li* (ki =9s,1,2,3,4,5), 
d+°Li* (kı = gs,1), 

3He + a, 


p+° Li’ Be* > 


(20) 
where kı denotes the energy levels of the first residual nu- 
clei, Mı, and gs denotes the ground state. Their energy-level 
schemes are taken from the experiments in Refs. [16, 45, 46]. 

For the first particle emission channel °Li(p, n)°Be*, the 
first residual nucleus “Be* that reaches the first energy level 
can still emit a proton with the second residual nucleus °Li*. 


Furthermore, the second residual nucleus °Li* can emit a pro- 
ton and alpha through the direct two-body breakup process, 
thus contributing to the (p, n2pa) reaction channel. 

For the first particle emission channel °Li(p, p)°Li*, the 
second excited energy level (Ek;=2 = 3.563 MeV) of the 
residual nucleus °Li cannot emit any particle because the par- 
ity is not conserved; thus, this reaction process purely con- 
tributes to the inelastic scattering channel. The first and third 
excited energy levels of Li can emit a deuteron, so they con- 
tribute to the (p, pda) reaction channel. If the first residual 
nucleus ®Li* is in the k-th (kı > 4) excited energy level, 
some energy levels will emit a proton with the second resid- 
ual nucleus °He*. Furthermore, the second residual nucleus 
5He* can emit a neutron and alpha from the direct two-body 
breakup process, thereby contributing to the (p, n2pa) reac- 
tion channel. If the first residual nucleus °Li* is in the k-th 
(kı > 4) excited energy level, some energy levels can also 
emit a deuteron, so these reaction processes contribute to the 
(p, pda) reaction channel. Therefore, the first particle emis- 
sion channel °Li(p, p)°Li* can contribute to the (p, n2pa) and 
(p, pda) reaction channels in the final state besides the elastic 
and inelastic channels. 

For the first particle emission channel, °Li(p, d)°Li*, the 
reaction process °Li* — p + a occurs as mentioned above, so 
this reaction channel belongs to the (p, dpa) reaction channel 


pound nucleus to the ground state of the first residual nucleus 
a, this process only contributes to the (p, 3He) a reaction 
channel. 


According to the analyses of the reaction channels dis- 
cussed above, the total spectra can be determined by adding 
all the partial spectra of the same outgoing particle obtained 
from every reaction channel. The contributions to the double- 
differential cross-sections of the total emitted protons are 
from elastic scattering, inelastic scattering, direct three-body 
breakup, and the (p,n2pa) and (p,pdq) reaction channels. The 
contributions to the double-differential cross-sections of the 
total emitted deuterons are from (p, d)°Li*, (p, pda), and 
the direct three-body breakup process ’Be* — p+ d + a. 
The contributions to the double-differential cross-sections of 
the total emitted alpha are only from (p,*He)a, (p, n2pa), 
(p, pda), and the direct three-body breakup process “Be* > 
p+d +a. The contribution to the double-differential cross- 
section of the total emitted *He is only from the (p, ?He)a 
reaction channel. It is worth mentioning that the experimen- 
tal double-differential cross-sections are a mixture of 7He and 
a because of the difficulties encountered in the measurement 
[22]. 


In conclusion, for the proton-induced 6Li reaction, reaction 
channels finally exist at an incident energy E, < 20 MeV as 


in the final state. For the first outgoing He from the com- follows: 
| 
onas J b= 95 (p, n) 
n+ Be ta p+°Li* +p+a_ (p,n2pa) 
ky = gs Compound elastic 
. 4+ 6Li* ky =2 (p, p’) 
p +îLi > Be —> 4 P ky =4,5 p+řHe*—>n+a (p,n2pa) (21) 
ky = 1,3,4,5 d+a (p, pda) 
d+°Li* k, = gs,1 pta (p, pda) 
3He +a (p, He) 
p+d+a three-body breakup (p, pda) 


IV. CALCULATED RESULTS 


The experimental double-differential cross-sections of 
light charged particles (proton, deuteron, “He, and a) for the 
p+°Li reaction were measured in Refs. [22, 23]. Based on 
the STLN model, a PLUNF code for the pt°Li reaction is 
developed to calculate the double-differential cross-sections 
of the outgoing nucleon and light composite charged parti- 
cles. Comparisons are performed between the calculations 
and measurements of the double-differential cross-sections 
for the total outgoing proton, deuteron, He, and alpha par- 
ticles for the p+°Li reaction. 


Because the target is contaminated by 'H, ‘Li, !°C, 
and t60, there are additional contributions to the double- 
differential cross-sections. Using Eq. (18), the kinetic en- 
ergies of the outgoing proton and its residual nuclei from 
the ground state to the fifth energy level of the pure Li and 


contaminants are identified. For example, for the double- 
differential cross-sections of the outgoing proton at an angle 
of 60° at Ep = 14 MeV in the LS, as shown in Fig. 2 (a), 
the black narrow bands in the below panel (b) represent the 
contributions of the outgoing proton from the discrete levels 
(only illustrated from ground state to the fifth energy level, as 
marked) of targets 6Li and contaminants (including 1H, 7Li, 
12C, and 160), From Fig. 2, one can see that the contribu- 
tions to the first abrupt peak on the right-hand side are the 
elastic scattering of protons with 1°C and 160. The contribu- 
tions to the second peak on the right-hand side are the elastic 
scattering of protons with °Li and ‘Li, and the first inelastic 
scattering of Li. The contribution to the third peak on the 
right-hand side is the pure inelastic scattering of protons with 
the first excited energy level of °Li. The contributions to the 
first peak on the left-hand side are the elastic scattering of pro- 
tons with 1H, the inelastic scattering of protons from the third 
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Fig. 8. (Color online) Same as Fig. 7, but for different outgoing 
angles, which are marked in the figure. 


to the fifth excited energy level of !?C, and the fifth excited 
energy level of “Li. The levels of “Li exhibit clear differences 
in their contributions to the double-differential cross-section 
of the outgoing protons compared with those of ®Li. 
Comparisons of the calculated double-differential cross- 
sections of the total outgoing proton with the measured data 
are shown in Figs. 3-5 at an incident proton energy of 14 
MeV for outgoing angles of 20°, 30°, 40°, 50°, 60°, 70°, 
80°, 90°, 100°, 110°, 120°, 130°, 140°, 150°, 160°, and 
165°. The black points represent the experimental data ob- 
tained from Ref. [23], and the red solid lines denote the 
calculated total double-differential cross-sections of the out- 
going proton. One can see that the calculated results agree 
well with the measurements, except for some peaks contam- 
inated by elastic and inelastic scattering from 1H, Li, 12C, 
and '°O, as reported in Ref. [23]. Taking the calculated 
double-differential cross-sections of the outgoing proton with 
60° at Ep = 14 MeV as an example, the partial spectra are 
shown in Fig. 6. The yellow dash-dotted line denotes the 
partial spectra of the proton emitted from a direct three-body 
breakup through ’Be* — p+d+a. The pink dash-dotted lines 
denote the partial spectra of the first outgoing proton from the 
compound nucleus to the ground state and the fifth excited en- 
ergy levels of the first residual nucleus, Li”. The green dash- 
dotted line represents the contribution of the reaction channel 
(p, np)°Li* — (p,n2p)a from the first excited energy level 
of °Be* to the first excited energy level of °Li*, which can be 
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Fig. 9. (Color online) Same as Fig. 3, but for the outgoing deuteron. 
The experimental data are taken from Ref. [22]. 


broken up into d+a. The brown dash-dotted lines denote the 
partial spectra of the second outgoing proton from the fourth 
and fifth excited energy levels of ĉLi* to the ground state of 
5He*. The purple dash-dotted lines denote the partial spectra 
of the second outgoing proton from the ground state and the 
first excited energy levels of *Li* to the ground state of +He. 
As the target is contaminated by 1H, ‘Li, 12C, and 160, there 
are additional contributions to the double-differential cross- 
sections besides the target nucleus Li. The energy levels 
of the contaminants contribute more to the lower outgoing 
energy regions as the outgoing angle increases, resulting in 
theoretical underestimations of the outgoing proton in these 
regions. 


The calculated  double-differential cross-sections 
of the total outgoing deuteron for the p+®°Li reac- 
tion at 14 MeV are compared with the experimen- 
tal data obtained at angles of 20°,30°, 40°, 50°, 60°, 
70°, 80°, 90°, 100°, 110°, 120°, 130°, 140°, 150°, 160°, and 
165°, as shown in Figs. 9-11. The black points represent 
the experimental data obtained from Ref. [22], and the red 
solid lines denote the calculated total double-differential 
cross-sections of the outgoing deuteron. Figure 12 shows 
the partial double-differential cross-sections of the outgoing 
deuteron with an angle of 60° at Æp = 14 MeV in the 
LS. The yellow dash-dotted line denotes the partial spectra 
of the emitted deuteron from a direct three-body breakup 
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Fig. 10. (Color online) Same as Fig. 9, but for different outgoing 
angles, which are marked in the figure. 
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Fig. 11. (Color online) Same as Fig. 9, but for different outgoing 
angles, which are marked in the figure. 


through “Be* —> p+d+a. The green dash-dotted lines denote 
the partial spectra of the first outgoing deuteron from the 
compound nucleus to the ground state and the first excited 
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Fig. 12. (Color online) Same as Fig. 6, but for the outgoing deuteron. 
The yellow dash-dotted line denotes the partial spectra of the emitted 
deuteron from direct three-body breakup through ’Be* — p+d + a. 
The green dash-dotted lines denote the partial spectra of the first out- 
going deuteron from the compound nucleus to the ground state and 
the first excited energy level of the first residual nucleus, °Li*. The 
pink dash-dotted lines denote the second outgoing deuteron from the 
first and third—fifth excited energy levels of °Li* to the ground state 
of “He. 


energy levels of the first residual nucleus, °Li*. The pink 
dash-dotted lines denote the second outgoing deuteron from 
the first and the third—fifth excited energy levels of Li” to 
the ground state of +He. 


The calculated double-differential cross-sections of the to- 
tal outgoing ĉHe and a for the p+®Li reaction at 14 MeV are 
compared with the experimental data obtained at angles of 
10°, 20°, 30°, 40°, 50°, 60°, 70°, 80°, 90°, 100°, 110°, and 
120°, as shown in Figs. 7—8. The black points represent the 
experimental data obtained from Ref. [22], and the red solid 
lines denote the calculated total double-differential cross- 
sections. The calculated results agree well with the measure- 
ments. From Fig. 13, one can see that the yellow dash-dotted 
line denotes the partial spectra of the emitted a from the direct 
three-body breakup process through “Be* —> p+d+a. The 
green dash-dotted line denotes the contribution of the reac- 
tion channel (p, np)°Li* — (p, n2p)a from the first excited 
energy level of ®Be* to the first excited energy level of °Li*, 
which can be broken up into d+a. The pink dash-dotted lines 
denote the second outgoing a from the first and third—fifth 
excited energy levels of ®Li* to the ground state of *He. The 
blue dash-dotted lines denote the partial spectra of the first 
outgoing He from the compound nucleus to the ground state 
of the first residual nucleus a. One can see that the first and 
second peaks on the right-hand side are the differential cross- 
sections of the outgoing *He and a for p+ °Li, respectively. 
In particular, p+°Li —++°He + a is different from a secondary 
particle emission because it has no threshold. It is worth men- 
tioning that the peaks are marked incorrectly in Figs. 7-8 
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Fig. 13. (Color online) Same as Fig. 6, but for outgoing “He and 
a. The yellow dash-dotted line denotes the partial spectra of the 
emitted a from direct three-body breakup through ‘Be* > p+d +a. 
The green dash-dotted line denotes the contribution of the reaction 
channel (p, np)°Li* — (p, n2p)a from the first excited energy level 
of °Be* to the first excited energy level of °Li*, which can break up 
into d+a. The pink dash-dotted lines denote the second outgoing a 
from the 1st and 3rd—Sth excited energy levels of °Li* to the ground 
state of “He. The blue dash-dotted lines denote the partial spectra 
of the first outgoing He from the compound nucleus to the ground 
state of the first residual nucleus a. The purple dash-dotted lines 
denote the partial spectra of the second outgoing a from the ground 
state and first excited energy levels of ŽLi* to the ground state of 
He. 


of Ref [22]. The purple dash-dotted lines denote the partial 
spectra of the second outgoing a from the ground state and 
the first excited energy levels of °Li* to the ground state of 
tHe. 


As an example, Fig. 14 shows the predicted total and 
partial spectra of the outgoing neutron for the p+ĉLi reac- 
tion at 14 MeV at an angle of 60°. The red line shows the 
total double-differential cross-sections of the outgoing neu- 
tron. The green dash-dotted lines denote the partial spectra of 
the first outgoing neutron from the compound nucleus to the 
ground state and the first excited energy level of the first resid- 
ual nucleus, ĉBe*. The brown dash-dotted lines denote the 
contribution of the reaction (p, 2p)°He — (p, n2p)a from the 
fourth and fifth excited energy levels of °He* to the ground 
state of a. 


However, the threshold energies of composite charged par- 
ticles, such as deuteron and triton (except He and alpha for 
the p+°Li reaction), are usually higher than those of the pro- 
tons for both the target and contaminants. Furthermore, the 
Coulomb barriers of these outgoing composite charged parti- 
cles are higher than those of the protons; therefore, the impact 
of the contaminants is greatly reduced. These factors hinder 
the accurate measurement of these composite charged parti- 
cles, particularly at larger outgoing angles. Thus, it is reason- 
able to assume that the theoretical calculations are slightly 
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Fig. 14. (Color online) The predicted total and partial double- 
differential cross-sections of the outgoing neutron from the p+°Li 
reaction with an outgoing angle of 60° at Ep = 14 MeV in the LS. 
The green dash-dotted lines denote the partial spectra of the first out- 
going neutron from the compound nucleus to the ground state and 
the first excited energy level of the first residual nucleus, °Be*. The 
brown dash-dotted lines denote the contribution of the reaction (p, 
2p)°He — (p, n2p) a from the fourth and fifth excited energy levels 
of °He* to the ground state of a. 


higher than those of the measurements for large outgoing an- 
gles. 


V. SUMMARY 


The STLN model [39, 47, 48] is improved to calculate the 
double-differential cross-sections of outgoing neutrons, pro- 
tons, deuterons, ?He, and alpha particles for the p+°Li reac- 
tion. The emission processes between discrete energy levels 
with energy, angular momentum, and parity conservations are 
strictly considered. The results show that the pre-equilibrium 
emission process is the dominant reaction mechanism for 1p- 
shell light nucleus reactions and that the calculated double- 
differential cross-sections are sensitive to the energy, spin, 
and parity of the discrete levels for both the target and resid- 
ual nuclei. Recoiling effects are also considered owing to the 
light mass of the 1p-shell nuclei. A PLUNF code is developed 
based on the STLN model to obtain an ENDF-6-formatted file 
of the double-differential cross-sections of the nucleon and 
light composite charged particles for the p+°Li reaction. The 
calculated results agree well with existing measurements of 
outgoing protons, deuterons, He, and alpha particles. 
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